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Abstract 

An equation previously proposed to describe the evolution of vortex line density in 
rotating counterflow turbulent tangles in superfluid helium is generalized to incorporate 
nonvanishing barycentric velocity and velocity gradients. Our generalization is compared 
with an analogous approach proposed by Lipniacki, and with experimental results by Swan- 
son et al. in rotating counterflow, and it is used to evaluate the vortex density in plane 
Couette and Poiseuille flows of superfluid helium. 

1 Introduction 

Many researches of quantum vortices in superfluids have been carried out on rotating systems 
and counterflow situations, both of them with vanishing barycentric velocity gradient pQ-[3]. 
Evolution equations have been proposed to describe the influence of heat flux and of angular 
velocity on the vortex dynamics [4] generalizing the well-known Vinen's equation for non- 
rotating systems |Tj — [3j [5]. An interesting challenge is to generalize these vortex evolution 
equations to include the influence of barycentric flow, which has much practical interest, for 
instance, in cryogenic applications. Here, we carry out such a generalization and we examine 
a recent proposal by Lipniacki [6j, which opens an interesting perspective but which, on the 
other side, discloses some aspects which have not been yet settled out with enough clarity. 

The aim of this paper is to generalize a previous equation proposed for rotating counterflow 
superfluid turbulence [4] by emphasizing more explicitly the dynamical role of the rotational 
of the superfluid velocity v s , related to quantized vortices. This allows us to write a proposal 
for the evolution equations of vortices in plane Couette and Poiseuille flows. In Section [3] 
we review some aspects of rotating counterflow and compare our generalized expression with 
Lipniacki's proposal [6|, which underlines the role of the polarization rather than of rot v s 
itself, and we stress some open problems. In Section H] we use a thermodynamic formalism to 
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relate the dynamical equation for vortices with a new term appearing in the mutual friction 
force, which we use for a comparison with that one by Lipniacki. In Section [5] we discuss 
several aspects of Couette and Poiseuille flows of superfluid helium including the presence of 
quantized vortices. 



2 Rotational of superfluid velocity and the dynamics of vortex 
line density 

An evolution equation for the dynamics of quantum vortices in rotating helium under coun- 
terflow was proposed in [4], describing the influence of the heat flow and of angular velocity 
on the vortex line density. In particular, the vortex-line density L was assumed to obey the 
following equation 



dL 



-(3kL 2 + 



L 3/2 



L. 



(2.1) 



where (3, a\, fo, Pi, and (3$ are dimensionless coefficients, k = h/m is the quantum of vorticity 
(to the mass of the 4 He atom and h Planck's constant), V = |V| (with V = v n — v s ) is the 
counterflow velocity, the relative velocity between averaged normal and superfluid velocities, 
which is proportional to the heat flux across the system, and f2 = |fi| is the angular velocity 
of the container. The values of the coefficients were obtained in [3] by comparison with 
experimental data of [7] and they were seen to satisfy the relations /?4 = \f2a.\ and (3\ = 
v2/32 — 2/3, which are required on relatively general arguments about the form of solutions. 
The values of the coefficients appearing in (|2.ip were independently calculated in [8j, and 
agree with those obtained in [3]. When f2 = 0, equation (|2,ip reduces to the well-known 
Vinen's equation [5], with parameters a\ and being respectively related to the production 
and destruction of vortices per unit volume and time. 

In [2J [9] it was shown that the value of coefficient a\ depends on the angle between the 
counterflow velocity V and Schwarz's binormal vector I [10] (see equation (|3.6p ). As observed 
in [2], ai = ayl ■ V, with ay the coefficient appearing in Vinen's equation (pure counterflow) 
PQ. Schwarz derived Vinen's equation using the vortex filament model obtaining ay = ac\, 
where a is the well-known coefficient appearing in the expression of the mutual friction force 
between vortex lines and the normal fluid and c\ denotes the average curvature of the tangle 
(see equation (13.7P ), In [4] in the regime of high rotation, the value I • V = 1/2 was found, 
so indicating that the vortex tangle is highly polarized. Coefficient f3 is linked to the average 
squared curvature of the vortices as (3k = a(3c^, with c 2 , defined in equation (13. 7p and (3 the 
vortex tension parameter, defined as ey = np s (3 with ey the energy per unit length of vortex 
line pQ. 

These equations lack an important source of vorticity, namely a barycentric velocity gra- 
dient, which is known to produce turbulence in many actual flows. Thus, it would be useful 
to generalize (|2.ip by incorporating in it barycentric velocity gradients. A possible way to do 
so would be simply adding new terms basing on dimensional analysis and on comparison with 
the observed phenomenology. Instead of proceeding in this way, we will interpret (|2.ip in some 
deeper terms, which will be useful for a consistent incorporation of the velocity gradient. 

To generalize equation (|2.ip we note that in the particular case of pure rotation £1 is related 
to rot v s as 2VL = |rot v s |, where v s is the macroscopic superfluid velocity. As Lipniacki noted 
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in a different proposal [6], writing an equation such as (J2JJ) in terms of rot v<j and V rather 
than in terms of f2 and V would be more general, because it would reduce to (12. ip for rotation, 
and it could be applied to other flows as plane Couette or Poiseuille flows (see Section [5]), where 
| rot v s | = dv sx (z)/dz, x being the direction of the fluid motion, z the direction orthogonal to 
the parallel plates, and v sx (z) the macroscopical superfluid velocity, depending only on z. 
In this way, the natural generalization of (|2.ip would be to rewrite it in terms of rot v s as 



dL n t2 

Tt = + 



<%\V + -^%\/ /drot v s 
V2 y 



L 3/2 



— rot v, 
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V2 



rot v« 



L. 



(2.2) 



Equation (|2.2p reduces to (|2.ip for pure rotation. Besides that, expression (|2.2p generalizes 
(|2.ip also on dynamical grounds. Note, indeed, that in (|2.1|) it is assumed that |rot v s | is 
equal to 20,. However, it will take some time for v s to get these values, by starting after some 
arbitrary initial state. Thus, whereas f2 in (|2.ip is taken as an externally fixed parameter, in 
(|2.2p rot v s is a dynamical quantity, which must be described by a suitable evolution equation. 
Then, the form (|2.ip will be useful after some transient interval, whereas (|2.2p is expected to 
be valid also for fast changes in v s . Further, equation (|2.2p can be applied also in different 
situations, as plane Couette and Poiseuille flows. Thus, equation (j2.2[) is the central point of 
this paper, as it generalizes (|2.ip both to a wider set of external conditions and to a wider 
domain of dynamical variations. 

Comparison with a similar approach by Lipniacki [6] will be useful for a better understand- 
ing of both approaches. Lipniacki [6j has essentially proposed to use as variable the so-called 
"polarity vector" (see also ), an important quantity in vortex dynamics, which he linked 
to the rotational of the averaged superfluid velocity 



P 



< s' > 



/ s'd£ _ V x v s 
fd£ ~ kL 



(2.3) 



For pure rotation one has p = Ct and V x v s = 2fi; therefore, we may rewrite equation 
(|2.ip in terms of p. Note that |p| £ [0, 1] measures the directional anisotropy of the tangent to 
the vortex lines: in particular |p| = 1 for a system of parallel vortices and |p| = for isotropic 
tangles. Thus, it is possible to express (|2.2|) in terms of p and to group the terms in it in a 
slightly different way, namely, in two groups, one of them with the factor VL 3 / 2 and the other 
one with kL 2 , mimicking in some way the form of the original Vinen's equation. In this way, 
we rewrite (12.21) as 



(3 V kL (3 kL 



+ aiVL 3/2 




(2.4) 



which, recalling (3\ = V202 — 2(3 and the previously mentioned relation 2Q = |rot v s | implying 
20/ kL = |p|, (|2.4p assumes the more compact form 



dt 



PkL 2 
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IPI 



(2.5) 



where B 



§h and A 



Pa 



In HI, coefficient B was found to be 0.89 while coefficient A is 



not properly a constant but undergoes a small step from 1 to 1.004 at the first counterflow 
critical velocity V c \. In this work, as already pointed out, we neglect this step assuming A = 1. 
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When inhomogeneities in the line density L are taken into account, the evolution equation 
for line density L must include a vortex density flux J L [12J 

§ + V-J L = a L , (2.6) 

where ctl stands for the production term given by the right-hand side of equation (|2.5p . The 
form of 3 L contains a connective contribution, Lvl with vl the velocity of vortex lines with 
respect to the laboratory frame, and a diffusive contribution. In some situations, when the 
rate of variation of the perturbations is higher than the reciprocal of the relaxation time of 
the diffusive flux [121 [13], one must take J L as an independent variable [H]. Here, neglecting 
the relaxation time of 3 L and considering isothermal situations, we take for J L the following 
simple law, where the diffusive contribution is analogous to Fick's diffusion law 

J L = -DVL + Lv-l. (2.7) 

The coefficient D (of the order of n p2],[T3]) is the diffusion coefficient of vortex lines. 

For a general hydrodynamic description, the evolution equations for v n and v s are needed. 
In particular, the evolution of v s is necessary to describe the evolution of rot v s in equation 
(|2.2|) . A set of equations frequently used are the Hall-Vinen-Bekarevich- Khalatnikov equations 
[TJ [15] , which in an inertial frame are written as 

Pr^T + P"(v„ • V)v n = - P —Vp n - p s SVT + F ns + ??V 2 v n , (2.8) 
at p 

Ps^ + Ps(v s ■ V)v s = -^V Ps + p s SVT - F ns + Ps T. (2.9) 
at p 

Here, p n and p s are effective pressures, defined as Vp n = Vp + (p s /2)W 2 , Vp s = Vp - 
(p n /2)VV 2 , p the total pressure, S the entropy, r\ the dynamic viscosity of the normal compo- 
nent, and p s T the vortex tension force, which vanishes for rectilinear vortices and for isotropic 
vortex tangles, but which may be relevant in other situations. In the situations considered in 
this paper, we will have T = 0. 

To describe the motion we need an expression for F ns , the mutual force between normal 
and superfluid components. The usual expression by Hall, Vinen, Bekarevich and Khalatnikov 

is ra 



ap s KL 



p x [p x (V - v;)] + — p x (V - Vi ) 

a 



(2.10) 



with q and a' being friction coefficient depending on temperature, and Vi the "self-induced 
velocity", which in the HVBK equations is approximated by 

v; = pv x p. (2.11) 

The expression for F ns must be consistent with the dynamics of L. In Section[4]we will explore 
how (|2.10p should be modified in order to be consistent with the evolution equation (|2.5h . and 
in Section [5] we will combine the equations in an analysis of plane Couette and Poiseuille flows 
in steady conditions. 
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3 Rotating counter flow 



In this Section, we investigate the proposed equation (12. 2ft for a rotating superfluid helium 
inside a cylindric container in the absence and in presence of counterfiow and we compare 
some results of our proposal with those of Lipniacki [6], and with the experimental data of 
Swanson et al. [7j. These authors considered a rotating container filled of helium II with an 
external counterfiow V parallel to the angular velocity ft of the container. For high angular 
velocities, they observed two critical counterfiow velocities V c \ and V c such that for < V < V c 
the line density L is approximately independent of V, undergoing only a small step (about 
0.4%) at the first critical velocity V c \ whereas for V > V c the line density L grows with V . 
Here we will neglect the small variation of L at the first critical velocity V c \, because our 
proposal reduces to (|2.ip in this situation — which was already carried out in [3] — and it is 
not necessary for the comparison with Lipniacki's proposal because the latter is valid only for 
V>V C . 

3.1 Pure rotation 

First of all, we consider the simplest situation of a cylindric container rotating around its axis. 
It is known that when the angular velocity O exceeds a critical value O c and the stationary 
state is reached, vortex lines parallel to the rotation axis are present whose number density 
follows the law L = 20/ k. The presence of these vortices may be explained observing that when 
the container begins to rotate the viscous normal fluid rotates with it, whereas the superfluid 
remains initially at rest, due to its vanishing viscosity. In this situation, the difference between 
v n and v s is zero along the rotation axis, but it is maximum near the walls of the container, 
that is, the counterfiow velocity increases for increasing distance from the axis. In this way 
the remnant vortices, which are formed during the cooling of helium and which are pinned to 
the walls, are influenced by the counterfiow velocity. This implies the growth of these vortices 
in agreement with the dynamical description proposed by Schwarz. According to this idea, 
vortices will grow near the walls, due to the relative velocity between normal and superfluid 
velocity, and will migrate towards the bulk of the system, forming in the stationary situation 
a regular array of vortices parallel to the rotation axis. 

The presence of vortices couples the normal fluid and the superfluid through the mutual 
friction force so that vortices are dragged by the normal fluid, and the average superfluid 
velocity v s becomes different from zero. This fact justifies the relation V x = 2f2 and the 
substitution of 20/kL = |p| in equation (|2.ip . At the light of the new arguments, in the case 
of pure rotation the vortex line density becomes L = 20/ k, which implies |p| = 1. 

Consider now equation (|2,5p in the case of pure rotation, when the stationary solution is 
reached, that is V =< |v n — v s | >~ 0. In this case (|2.5p has two stationary solutions, |p| = 1 
and |p| = 1/B 2 . As one can easily verify, the solution |p| = 1 is stable if B < 1 and this is 
the case because the coefficient B was found to be 0.89 [4j. To describe the non-stationary 
regime, one needs to introduce equations (|2.8p and (12.91) for the averaged normal and superfluid 
velocities. 

3.2 Fast rotation and external counterfiow 

Equation (12. 5|) can be also written as 
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dL L 3 / 2 { 1 - ^\p~\) \ ai V - ^L 1 ' 2 ( 1 - By^\)] . (3.1) 



dt 

As we have pointed out above, pure rotation is well described by (|3.ip . because in this situation 
v s = and |p| = 1 is a stationary solution of (|3.ip . meaning complete polarization. The 
non-zero stationary solutions of (|3.ip are 



|p| = l and L^ = ^V + B\ ] ^^. (3.2) 

\3k V k 

To study the stability of the solution |p| = 1, we linearize Eq. (|3.ip for the perturbations. In 
the hypothesis that the perturbation 5 does not modify the vorticity to = rot v s , the relation 
<5|p| = — (\p\/L)5L is obtained, which allows us to obtain the following evolution equation for 
the perturbation SL 



dSL \ 



n|1 \pk{i-b)l 



2LV2 2' 



6L. (3.3) 



IpI=i 

From the previous equation it follows that the solution |p| = 1 is stable for V less than 

Vc = — (1-S)V|Vxv s |k, (3.4) 

which corresponds to the critical velocity V c in the experiments of Swanson et al. [7]. Note 
that if B = 1 in Eq. (|3.4I) . the critical counterflow velocity for which the straight vortex lines 
parallel to the rotation axis become unstable is zero. From an experimental point of view 
this is not the case because a nonvanishing critical velocity is observed, confirming the value, 
B = 0.89 < 1, obtained in reference [JJ. 

For counterflow velocity higher than the critical velocity (|3.4p . the solution |p| = 1 becomes 
unstable, and the line density L assumes the value p. 2b ) which depends on V and |rot v s |. 

Now, we consider the second term in the right hand side of (|3.2b ). namely B\l i^iXfl 
For low values of the counterflow velocity, the vorticity is essentially due to the rotation, and 
therefore we put |V x v s | = 2fi, recovering the results obtained in [JJ. 

3.3 Comparison with Lipniacki's proposal 

Recently a hydrodynamical model of superfluid turbulence was proposed by Lipniacki [6], 
mainly with the aim to studying the hydrodynamics of partially polarized tangles arising in 
rotating counterflow or in plane Couette flow. Thus, it is interesting to compare with his work, 
whose aims are similar to ours. 

Lipniacki writes Vinen's equation as 

^ = aL 3 / 2 Cl (|p|)I • V - (3a 2 c 2 (\p\)L 2 , (3.5) 

where j3 is a constant of the order of k, and a the friction coefficient appearing in the expression 
of the mutual friction force; I is the binormal vector, 

< s ; x s" > 

1 = <|s"|> ' (3 - 6) 



(3 



defined by Schwarz [10] to describe the polarization of the binormal s' x s", of the vortex lines, 
with s' and s" being the first and second derivatives of the curve s(£) describing a vortex line 
with respect to the arc-length £, s' the unit tangent along the line and s" the curvature vector. 

The coefficients c\ and c 2 . measure the average curvature and curvature squared of the 
tangle, respectively. They are given, according to the microscopic model by Schwarz [10J, by 

s \d£, 



ci 



4 



AL 2 



J'\2 



(3.7) 

Lipniacki proposes 

that ci and e 2 . should depend on the polarization |p| , and that they should vanish for completely 



AL 3 /2 

where A is the volume on which one makes the averaging indicated in ([3 

-2 



polarized tangles because in this case s" = for all the vortex lines. To describe the reduction 
in ci and c 2 , with respect to its usual variable for a nonpolarized tangle, which will be designed 
as cio and c| , respectively, he assumes that 

|21 2 



ci(|p|) - cio [1 - |p| 2 ] , c|(|p|) ~ 4 [1 - |p| _ 
In contrast, our expression (|2.5p could be interpreted in this perspective as 



ci(|p|) - c io 



iPl 



cI(IpI) 



- C20 



IPl 



1 - 



IPl 



(3. 



(3.9) 



Therefore, it arises the question of the comparison of both equations (|2.5p and (|3.5p with the 
experimental data, and a deeper understanding of the influence of polarity on the coefficients 
ci and c|. 
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Figure 1: Comparison of the stationary solutions of Lipniacki's model (I3,10p (dashed line) 
and Jou and Mongiovi's model (|2.5p (black line) with the experimental data (solid circles) by 
Swanson et al. for counterffow velocity bigger than the second critical velocity V c and angular 
velocity 0.2 Hz, 0.4 Hz, 0.6 Hz, 0.8 Hz and 1 Hz. Lipniacki's model does not give the horizontal 
part of the plot, corresponding to V < V c . 

The evolution equation for the vortex line density L, proposed by Lipniacki in [6], has the 
explicit form 

^ = aI,c w VL^ [1 - |p| 2 ] - /3«4L 2 [1 - |p| 2 ] 2 , (3.10) 
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where Io = I • V, and the subscript stands for independence of Jo on V and L. The author 
chooses for To the same values found in pure counterfiow, in such a way to not consider the 
anisotropy of the vortex tangle, which is present owing of the high values of rotation considered 
in the experiments by Swanson et al. [7j. 

Equation (13. lQf) . as the author remarks, does not describe any of the two critical velocities, 
V c i or V c , of the experiments of Swanson et al. [7j. Lipniacki's aim is instead to describe the 
relation between angular velocity, counterfiow, and line-length density for polarized tangles 
above the second critical velocity V c . This implies the need of a comparison, in the uniform 
steady rotation and counterfiow, between (j2.5|) . (|3.10p and the experimental data of Swanson 
et al. 0. 

The stationary solutions of the equation (|3.10p are |p| = 1 (which however is unstable) and 

L = - Lh n2 , (3.11) 



where 



1 - (L u /L) 



( ciojo \ 2 i t [rotyaj _ 2tt 



L H = V 2 [ and L ^ = = — (3-12) 

are the steady state vortex-line density in pure counterfiow and in pure rotation, respectively. 

In Fig. [H we compare the results of equations (|2.5p and (|3.10p with the experimental data 
of the Fig. 2 of Swanson's experiments. It follows that (|2.5p (black line) describes better the 
experimental data (solid circle) than (|3. 10[) (dashed line), not only for V > V c , but it also 
yields the horizontal branch of the results for V < V c , which are not described by equation 
(|3.1Up . Comparison with experimental data shows that in the considered range of values of V 
and f2 equation (|2.5p fits better the experimental results. 

A reason for the difference between proposals (|2.5[) and (I3,10p could be related not to 
the evaluation of the integrals in (|3.7p but to a different microscopical interpretation of some 
terms in the evolution equation for L. Schwarz's derivation [10| is based on the dynamics of 
vortex breaking and reconnection, and its production and destruction terms tend to zero for 
completely polarized systems, as rightly pointed out by Lipniacki. However, the origin of the 
rotational terms in (|2.ip could be completely different. It is known that in rotating superfluid 
helium the vortices grow near the walls due to the rotation, and drift towards the center of 
the system, where they find a repulsion due to other vortices. These forces are different from 
zero even for completely polarized vortices, in contrast to the terms from ()3.7|) . It could then 
be that the vanishing of the terms in (|2.5p as 1 — had a different physical origin than 

the vanishing proposal by Lipniacki from a different model. These open questions stress the 
need of the inclusion of rotational effects in a more general version of Schwarz's derivation of 
Vinen's equation. 



4 Thermodynamic analysis of polarized superfluid turbulence 

In this Section, we will perform two modifications of the expression of the mutual friction 
force, as used in the HVBK model, which are necessary to incorporate the anisotropy of the 
vortex tangle and to insure the thermodynamic consistency of the evolution equation for L 
and for v s , according to the formalism of linear irreversible thermodynamics [9j [16]. Since 
(|2.5p differs from the usual Vinen's equation, it is logical to ask how these modifications will 
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change the form of F ns . For the sake of simplicity, we will neglect here the contribution of the 
self- induced velocity in (|2,10j) . 

First, we will take into account the anisotropy of the tangle introducing the tensor II = 
IP + IP, studied in [9], p], 

IT = - < U-sV >, n a = -— <W-s' >. (4.1) 
2 la 

In this equation s' is the unit vector tangent to the vortex lines, sV is the diadic product, U 
is the unit matrix, W is the Ricci third-order tensor and the angular brackets stands for the 
average over vortex lines in a given volume. The tensor II s describes the orientation of the 
tangents s' of the vortex lines, and the tensor II a — associated to an axial vector — describes 
the polarization; in other words, II a is related to the first-order moment of the orientational 
distribution function of s' and II s is related to second-order moment. As shown in Ref. 
using tensor II, the mutual friction force can be written 

F ns = -ap s KL^n ■ V. (4.2) 

If we suppose isotropy in the tangle, it results II s = U, II a = and one finds the usual 
expression 

F ns = ~a Ps KLV. (4.3) 

The tensor II in (|4.2p allows one to deal under a same formalism an array of parallel 
straight vortices as well as an isotropic tangle, and also the intermediate situations. 

Now, we follow the general lines of [9], [T7] to propose a modification to (|4.2|) with the aim 
to determine an evolution equation for v s consistent with (|2.5|) . According to the formalism of 
nonequilibrium thermodynamics one may obtain evolution equations for v s and L by writing 
dv s /dt and dL/dt in terms of their conjugate thermodynamic forces — p s \ and ey- The evo- 
lution equation (|2.9p for v s , neglecting inhomogeneous contributions of pressure, temperature 
and velocity, in an inertial frame, is written 

p s ^- = -F ns = apsKL^n. ■ V. (4.4) 



However, in the right-hand side of (|4.2p must be included additional contributions to make 
(|4.4p thermodynamically consistent with (|2.5p . 

In a way similar to that presented in [9], we write dv s /dt and dL/dt in matrix form in the 
system (I4.5p . In it, we write the equation for L in the form given in equation (]2.5p and by 
means of Onsager-Casimir reciprocity we obtain an additional contribution to the evolution 
equation for v s . The result is 

*\-J -A a "i n ±-i" 2 (i-vW)v )(-?,v\ 

tr -?i" ! i-v*i v -±l i-j® ) i-bj® H * ) (4 ' 5) 



The sign ambiguity present in that equation comes in a natural way from the Onsager-Casimir 
reciprocity relation. Indeed, in Feynman-Vinen view, L is a scalar quantity which does not 
change under time reversal, unlike the superfluid velocity v s which changes sign. According 
to Onsager-Casimir, this leads to antisymmetry of crossed coefficients thus leading to the + 
sign. In Schwarz view, L possesses vectorial properties and it would change on time reversal, 
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just like the superfluid velocity. This leads to the symmetry of the kinetic coefficients in the 
matrix in (|4.5h . i.e. to the - sign in the upper right-hand term. Below, we will directly take 
the minus sign, for the sake of a more direct comparison with the work by Lipniacki. 
Therefore the equation for dv s /dt becomes 

Ps^- = ap s KLhl • V - eyarL 1 ' 2 (l - ylpf) V. (4.6) 

The new term not contained in the evolution equation (|4,4p for v s is the coupling term between 
dv s /dt and ey in the matrix in (|4.5p . Note that this term depends on the direction but not 
on the modulus of V. This class of terms are called dry-friction terms. 

Observing that in the steady state (L, |rotv s | and V constant) the solutions of vortex line 
density equation (|2,5p can be written as 



L y 2= /JrotV^ for 0<v<v ^ (4?) 



L 1 ^ 2 = -tt-(V — V c ) + a/^, for V > V c , (4.8) 
and substituting them in (|4.6|) . we obtain the following expression for the coupling force 

V = 0, for V < V c , (4.9) 



1 coupl — €yQL\ 



L l/2 I l rotv - 



' coupl €yCX\ 



L i/2 / 1 rot v.. 



V = ey^(V-V c )V, for V > V c . (4.10) 



As a consequence, for V < V c the coupling force is absent (as in pure rotation) while, for 
V > V c , when the array of rectilinear vortex lines becomes a disordered tangle, the additional 
term (14. 9p appears. Indeed, in a almost-steady state (L and |rotv s | constant), for V < V c , 
equation (|4.4p would be valid, with L expressed by (|4.7|) . whereas, for V > V c it would become 

f= aK L 2 -n-V + CV] ^- { V-V c )V, (4.11) 

with L expressed by ()4.8p . Summarizing, in steady states for V < V c the dry- friction force 
is absent, while it appears for V > V c , when the array of rectilinear vortex lines becomes a 
disordered tangle. Thus V c indicates the threshold not only of the vortex line dynamics but 
also of the friction acting on the velocity v s itself; this seems logical, as both variables are 
mutually related. 

Summarizing, in this Section we have proposed to substitute the expression (|4.3p of the 
mutual friction force used in the HVBK model with 

F ns = -ap sK L^n ■ V - eyaiL 3 / 2 (l - y1p|) V (4.12) 



for the sake of thermodynamic consistency with (|2.5p . 

To complete the comparison between Lipniacki's and our model, we analyze in both models 
the expression of the mutual friction force, which in HVBK equation is expressed by (|2.10p . 
while in general terms it is expressed as 

F ns = ap s nL < s' x [s' x (V - v } )] > +a p s nL < s' x (V - v ; ) > . (4.13) 
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Lipniacki neglects the coefficient a' in Eq. (14.131) . and he approaches the quantity < s' x (s' x 
V) >~ [< s's' > -U]V = I V -V (where L, =< s'(s' • V) >) with 

< s' x (s' x V) >~ p x (p x V) - -(1 - |p| 2 )V, (4.14) 

3 

and the quantity < s' x (s' x Vi) >~ (3 < s' x s" >= (3ciL l ' 2 l with 

< s' x (s' x vi) >~ -/3/ ci (l - |p| 2 )L 1/2 V. (4.15) 

In explicit terms he uses 



F ns = a.Kp s L 



p(p • v) - ^ 2 - ± y L + 0Wi - IpI 2 )£ 1/2 v 



(4.16) 



So in the work of Lipniacki, the tensor |n s = < U — s's' > assumes the expression: 

2 -Il s ~ [U - pp] + |(1 - |p| 2 )U = 5 ~3 IP|2 U - pp. (4.17) 

Note that (|4.17p does not respect the relation trace[< U — s's' >] = 2, following from the 
normalized character of s', if |p| ^ 1. In fact it is 



trace 



5 - 2|p 



2 



-U-pp 



3|p| • (4.18) 



The last term in (|4.16p is a consequence of the drift of the tangle in the direction of the 
counternow, caused by its anisotropy, where I = IoV. This term is substituted in our model 
by the last term in (|4.6p . which we can rewrite, recalling that ey = p s K,(3 and ot\ = ac\oIo as 

-eyQiL 3/2 (l - ^/\p~\) V = -p s K(3ac w I L 3/2 (l - y/\p\) V. (4.19) 



' coupl ' 

As it is seen, this term differs from the one of Lipniacki, in the contribution due to the 
polarization of the tangle, which in our approach depends on 1 — \/|p| , and in Lipniacki's one 
on 1 — |p| 2 . We note also that, in this interpretation, we must choose the negative sign in the 
expression of this coupling term, in agreement with the microscopic derivation of the filament 
model by Schwarz. 

Lipniacki does not consider the tension T. For a comparison with our work, we must 
observe that in Lipniacki's model the quantity < s's' > is approximated by pp, and this 
approximation is correct only if most of the vortex lines in the volume have the same direction. 

In Ref. [TT] we have provided a microscopic paramagnetic analogy to relate p =< s' > 
with ft and V, in the case of simultaneous counternow and rotation, but we have not studied 
the statistic of the curvature vector s". In contrast, Lipniacki leaves open the value of p and 
makes some simple hypotheses about < |s"| > and < |s"| 2 > in his analysis of the possible 
influence of polarization in the Vinen's equation. 

A further difference between our model and that of Lipniacki refers to the form of the 
vortex flux for which he writes 



J L = Lv L = L 



v s + ap x V + f3al c w {\ - |p| 2 )L 1/2 V + /3aI k L 1/2 l (4.20) 



where the vector Ik represent the curvature of Co s lines. This last term is exactly zero if the 
vortex lines are closed, isotropic of straight, and otherwise it is expected to be small, except 
for the case when all the vortex lines are parallel to each other but bent. This is only the 
convective contribution, to which it should be added the diffusive contribution 3 L = —DVL. 
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5 Vortex-line density in steady plane flows 



In equation (12.2p (and equation (]2.5p ) we have rewritten previous equation (|2.ip for rotating 
counterflow turbulence in liquid helium in terms of |rot v s |. For pure rotation, |rot v s | = 2f2 
and we just have our original equation, but (|2.2p may be also used to describe situations with 
barycentric motion as plane Couette and Poiseuille flows (without external heat flux) between 
two parallel plates. Here we will consider two plates separated by a distance D, one at rest 
and the other one moving at velocity Vo (Couette flow), or plane Poiseuille flow, given by 
a longitudinal pressure gradient along the direction of two parallel quiescent walls. Here, we 
will deal with steady states and quasi-stationary states. We will assume that the flow of the 
normal component remains laminar. This requires that the Reynolds number DVq/t], with rj 
the viscosity of the normal component and Vo the characteristic velocity of flow, is sufficiently 
small. On the other side, in analogy with the rotating container, we assume that the velocity Vo 
is sufficiently high to neglect the "effects of the walls" [18J. The essential problem in both cases 
is to find the distribution of the superfluid velocity, the vortex density and the mutual friction 
force. We will focus our attention mainly to steady state situations, as simple illustration of 
the changes implied by the new equations (|2.5p and (|4.6|) . for L and v s . 

5.1 Plane Couette flow 

We assume two plane surfaces at z = and z = D such that the second one moves parallely to 
the first one at the velocity Vo, and that the relative velocity between normal and superfluid 
velocities has a profile V = (V x (z),0,0). The dynamics of vortex formation is similar to that 
in the rotating cylinder. When the upper plate starts suddenly moving with respect to the 
lower plate, the normal component starts moving under the action of the viscous force and 
non-slip condition, whereas the superfluid component will remain initially insensitive to the 
motion of the plate. Thus, a relative velocity (the counterflow velocity) V = v„ — v s will arise 
between the two components. This counterflow velocity V depends on the distance from the 
lower plate, in particular V is maximum for z = D (upper plate) and zero for z = (lower 
plate) . 

When the counterflow velocity reaches a critical value near the moving plane, the remnant 
vortices, always present in He II, begin to lengthen and reconnect to form other vortices, which 
diffuse towards the lower plate (at rest) forming, in the stationary situation, an array of vortices 
parallel to each other and orthogonal to the flow. Through the vortices, the normal and the 
superfluid components become coupled by the mutual friction force F ns , and the superfluid 
will tend to match its velocity with that of the normal fluid (V — ► 0); this will introduce a 
rot v s ^ in the superfluid, expressed by \dv s /dz\. After a sufficiently long time, it is expected 
that a steady shear flow will have formed, with v n = v s depending only on z and having the x 
direction and such that dv n /dz = dv s /dz = Vq/D, corresponding to the physical Newtonian 
linear profile, which follows from (|2.8p and (|2.9p with vanishing tension force T = 0, and f)4. 12[) 
in which F ns = for V = and |p| = 1. Then, it results |rot v s | = \q/D. 

Introduction of this value in (|2.2p would give the areal density of parallel and straight vortex 
lines, perpendicular to the flow. However, as it has been said in Section 2, the replacement of 
0, in terms of rot v s is deeper than a formal substitution because v s will not become related 
to the gradient of the barycentric velocity until a complex transient process has lapsed. Thus, 
the direct replacement of 2f2 in (|2.ip by dv sx /dz in shear flows, with v sx the ^-component 
of the macroscopic superfluid velocity, will be valid for steady states and for relatively slow 
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variations with respect to steady states. Otherwise, rot v s with its own nontrivial dynamics 
should be considered in (12, 2j) . The situation of Couette flow may be rather illustrative of these 
features. 

Then, the dynamics of L in this case is described by 



dL _ r9 



dv s 



dz 



L 3/2 



2 



dv s 



dz 



dv s 



dz 



L-V-J L , (5.1) 



where the coefficients should obey the relations indicated below (12, ip . and where the last term 
stands for the effects of the vortex flux for inhomogeneous systems. 

In the stationary situation V ~ and, according to (15. ip . there will be a completely 
polarized array of vortices, perpendicular to the velocity and to the velocity gradient, with 
uniform areal density given by 



i = I 

K 



<9v. 



dz 



kD' 



(5.2) 



In this view, the stationary character of L would require V to be zero, in such a way 
that normal fluid, superfluid and vortices would move at the same speed and without internal 
friction. However, equation (|5.ip has the intrinsic feature that for V less than a value V c the 
vortex line density does not depend on V and is given by (|5.2p . This critical relative velocity 
is, according to (|5.ip . 



Vr 



£ 




(32 


i 






dv s 






p. 




dz 




dz 



(5.3) 



with d ~ 3.7, according to the values of the coefficients used in (|2.ip to describe the value of 
V c in rotating counterflow velocity. 

This indicates that the ordered array of vortices formed in the Couette flow is stable until 
V <V C . This means that, as Vb grows, the regular array of rectilinear vortices, orthogonal to 
Vo, is still present and the velocities v n , v s and V have rectilinear profiles, but with slightly 
different slope. The value of V is maximum near the plane z = D. When the counterflow 
velocity V reaches the critical value V c the regular Couette array of vortices becomes unstable 
and a disordered tangle of vortex lines appears between the two plates in the zone. If one 
did not apply (|2.1|) . but only intuitive reasoning without the detailed quantitative analysis 
showing this critical velocity, one would expect that for V > will always be a disordered 
tangle of vortices. 



5.2 Plane Poiseuille flow 

Equation (15. ip may be applied to plane Poiseuille flow between two quiescent parallel walls 
at z = dbJD/2, driven by a longitudinal pressure gradient. The steady velocity profile for a 
Newtonian viscous fluid is parabolic, and has the form 



V x (z) 



Ap 
2^1 



Pi 

4 



1 



4z*_ 
D 2 



V77- 



1 



Az^ 
D 2 



(5.4) 

The pressure 



with the pressure gradient, 77 the viscosity and V max = {D 2 Ap) / (877Z) . .. , 

gradient acts on each component in the proportion established by the HVBK equations (|2.8I) 

(ESI). 
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(a) (b) 

Figure 2: Initial profile (a) and steady profile (b) of the superfluid (continuous line) and normal 
velocities (dashed line). 

Initially, the velocity profile of the normal component, submitted to viscous effects and to 
no-slip conditions on the walls, will be rather different from that of the superfluid component, 
which may slip freely along the walls (see Figured). As a useful simplification, one may ap- 
proximate the velocity profiles as a parabolic (Poiseuille profile) and a flat profile, respectively 
|19j . which equals one to each other in two points at the distance zq from the center of the 
plates. Then, one must search how these profiles will evolve under their mutual interaction 
due to the friction force, caused by the presence of the vortices. 

During the transient regime, vortices will be produced mainly in the regions where the 
relative velocity V is higher than a critical value V c — which may be also the central region — 
but they will be transferred to the matching region where v n = v s because of the second term 
in expression (I2.10P for the mutual friction force, which is a Magnus force yielding a vortex 
lateral drift velocity described by vuiateraldrift) = a ' s ' x V. The accumulation of vortices in 
the region where the two fluids have the same velocity will enlarge the width of the matching 
region (the profile of v s tends to the profile of v n ), until arriving at a situation where V will 
be lower than V c so that not more vortices will be produced. The steady profile will have the 
approximate form of Figure [2b, similar to that considered by Samuels (Fig. 7 of |20j), but in 
the matching region v n and v s are not exactly equal, in contrast with Couette flow or rotating 
cylinder, because there is need of a friction force to cancel out the term in the pressure gradient 
in the HVBK equations, as shown in (|5.6f) below. 

In the steady state, for isothermal flow, and neglecting the tension T, which vanishes for 
rectilinear vortices and for isotropic tangles, equations (|2.8p - (|2.9p reduce to 

- — Vp n + F ns + r?V 2 v n = 0, (5.5) 
P 

-^Vp s -F ns = 0. (5.6) 
P 

By adding these equations one obtains —Vp + ryV 2 v n = 0, which shows that the velocity 
profile of the normal component is the usual one corresponding to the motion it would have 
by itself, without the interaction with the superfluid unless some contributions with T ^ 
would appear, in the form, for instance, of local anisotropy vortex tangles. On the other side, 
from (|5.6|) it is seen that F ns will be different from zero, given by F ns = — ^\/p s . Thus, v n 
and v s will be slightly different, if Vp is low enough, and there will be an array of straight 
vortices, which we calculate below. 
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The most relevant features of the steady profile are: the width 2z c of the central zone 
without vortices and flat v s profile, the width z w of the boundary layer also without vortices 
and flat v s profile, and L, the averaged vortex density in the matching zone where vortices 
concentrate. We will compute them from simple qualitative arguments. 

To compute z c and z w we will ask that the corresponding circulation of V ns in these regions 
is lower than the vorticity quantum k. This is a sufficient condition for the lack of vortices 
in this zone. The argument is similar to that which could be used to estimate the critical 
angular velocity for the formation of the first vortex line in a rotating cylinder. If the cylinder 
is rotating with angular speed the circulation of V will be 2ttR 2 Q,; to obtain Q c one equates 
this quantity to k, and one gets ft c = k/(2itR 2 ). The exact result is fl c = Kln(b/ao)/(2irR 2 ) 
PQ, with ao the radius of the vortex line and b a distance of the order of the line spacing, 
which in the case of the one vortex is of the order of the radius R of the cylinder. In the plane 
Poiseuille flow situation the quantity b is of the order z c , in the central zone, and of the order 
z w , in the boundary layer zone. 

Thus, to estimate z c we calculate the circulation of V ns = [z 2 — z 2 ~\ in the zone between 
z = and z = z c and equate it to Kln(b/ao). One has 



2 2" 

z c -z 



2 = 



Ap 

dx = YJ z c ~ Kln(cz c /ao) (5.7) 



where 7 is the contour of the square whose side is z c and c is a numerical constant of the order 
of the unity. This may be expressed in terms of the maximum velocity V max of v n as given by 
(|5.4|) . leading to expression 

zl n\n(cz c /a ) 



D 3 ADK 



(5.8) 



max 



Concerning the width of the boundary layer z w , a similar argument yields 
V ns -dl = 



71 




D , 2 



2 ' 



dx = — \Bz 2 w - z%] w Kln(c'z w /a ), 
D 2r ll 



1 



(5.9) 

where 71 is the contour of the square whose side is z w and d is a numerical constant of the 
order of the unity. Up to second order in z w , this yields 

— Dzl = K\n(c'z w /a ), (5.10) 
lf]l 

and using expression (|5.4p for the v n profile, the previous expression may be rewritten in terms 

of Vmax 3-S 

zl_ = Kln(c J z w /a ) 

D 2 4DV max ' 1 " ' 

This expression is similar to the one obtained by Samuels in [20] for the width of the outer 
layer in a cylindrical Poiseuille flow in a tube of diameter D (his eq. (15)), which was 

r c x2 



— ln(— ). (5.12) 

From (|5.8p and (|5.1ip it is found that the widths z c and z w decrease for increasing V max as 
z c ~ Vmax 3 and z w ~ Vmax' ■ Thus for increasing V max (i.e. increasing pressure gradient) the 
central zone and the outer zone boundary layer free of vortices will become thinner. The flat 
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profile of v s in these zones is consistent with the absence of vortices, according to the relation 
L = \dv s /dz\/ k, analogous to the expression (15. 2h . and which vanishes for flat profile. 

In the matching region the value of v n — v s will not be strictly zero, but because of 
restriction (15. 6j) if v n — v s is approximately constant in this region, one will have that 



Uz) 



1 


dv s 


~, 1 


dv n 


K 


dz 


K 


dz 



8V„. 



iD 2 



(5.13) 



It is known that there exist two values of z where the velocities, v n and v s , are equal, but, 
in general, in the rest of the z domain they could not be exactly equal. This implies that 
the mutual friction force could depend on z and that the distribution of the vortices could 
not be homogeneous. To overcome this problem, we average the value of L in the region 
between z = z c and z = z\ = D/2 — z w . Of course, the value of L in the region between 
z = —z\ = —D/2 + z w and z = — z c will be the same of the first region by symmetry. To 
estimate, we assume that the averaged profile of the superfluid velocity can be approximated 
by the profile of the normal velocity, so obtaining 



L 



8Vrr 



kD 



2V r , 



kD 



1 + 2 



Zc 

D 
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(5.14) 



Introducing z c and z w as obtained from (15. 8p and (|5.1ip we would have an estimate of L in 
terms of V max , or, equivalently, in terms of Ap. A more detailed analysis could be carried out 
from the transversal terms of the vortex flux, where the Magnus drift and the diffusion flux in 
(|2.7p would cancel each other. 

Expression (|5.8p may be used to obtain the conditions for a laminar flow without any 
vortex. This situation will be found when the width of the central zone without vortices z c is 
bigger than D/2. This leads to the condition DV max / 1 k < 2 In (D/(2ao)). For D ~ 10 _2 m, and 
since ao ~ lCT 10 m, we have DV max / n < 20. The dimensionless quantity DV max / 'k is analogous 
to the Reynolds number. In viscous fluid, the Reynolds number is defined as DV/v, with v 
being the kinematic viscosity u = rj/p, which has dimensions m 2 s _1 . The vorticity quantum k 
has also dimension m 2 s _1 and therefore DV max / ' n may be considered as a quantum Reynolds 
number. A similar number has been used in pure counterflow experiments in cylindrical 
containers of diameter D where, for instance, the appearance of the first vortex takes place 
at T = 1.7K for DV/k ~ 80 |21j. A more rigorous stability analysis would be convenient to 
obtain more values of the critical quantum Reynolds number in both situations. 



6 Conclusions 

The quantized character of vorticity in superfluids emphasizes the special importance of vortex 
lines, whose dynamics becomes a central aspect of rotating or turbulent flows of superfluids. 
The main proposal of this paper is equation (|2.2p for the evolution of vortex line density, which 
generalizes our previous proposal (|2.ip for rotating counterflow situations. Here, by writing the 
local average rotational of the superfluid component instead of the angular velocity, we have 
enlarged the set of applications of the theory in two main aspects. One of them is that (12. 2p . 
in contrast to (|2.ip . may be applied not only to rotation but also to shear flows, as illustrated 
in Section [5j The second enlargement is of dynamical nature: in (12. 2p rot v s appears, and 
v s itself should satisfy its own evolution equation, which is coupled to the evolution of v n , 
the velocity of the normal component. Then, (|2.2p becomes deeply coupled to the HVBK 
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equations (|2.8j) and (12. 9p for v n and v s not only through the mutual force F ns between the 
normal and the superfluid components, which requires the knowledge of L, but also because 
in (|2.2p v s is needed to obtain L. Thus, the coupling of these equations is much emphasized 
in (|2.2I) as compared to (|2,1D , 

For situations close to nonequilibrium steady states or for slow variations of v s , in such 
a way that rot v s is well described by the angular velocity or by the barycentric velocity 
gradient, equations (|2.ip or (|5.ip describe the vortex line density in terms of Q or dv sx /dz. In 
each case we have provided an estimation of the vortex density and of the superfluid velocity 
profile in the steady state. 

We have compared our proposal with that of Lipniacki, which shares the objectives of the 
present paper but stresses the polarization p = rot v s /kL more than rot v s itself. Lipniacki's 
evolution equation for L is, essentially, the classical Vinen's equation, but with the new aspect 
that its coefficients become dependent on the polarization p according to the microscopic 
identification of the coefficients proposed by Schwarz [10j . The disagreement between our 
equation (|2.5p and the Lipniacki's proposal (|3.10p may be due to the different physical origin 
of the terms dependent on the polarization. Our opinion is that Schwarz derivation of Vinen's 
equation (|3.5p does not include some relevant contributions of rotational systems. For straight 
parallel vortices, as those arising in pure rotation experiments, both the production and the 
destruction terms vanish. This is consistent with Schwarz's postulates for the vortices, but in 
purely rotational flows the dynamics of vortices has a different origin, related to the migration 
of vortices formed on the wall towards the center of the system, and with repulsion forces 
amongst vortices. Thus a general treatment would require to include these effects besides the 
Scharwz effects in (13.5p . and it could provide a further understanding of the differences between 
(|2.5p and (|3.10p . In any case, comparison with experimental results in Fig. 1 indicates that 
(|2.5p yields a better description of them. 
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